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Abstract 

We present a unified approach to constrained implicit Lagrangian and Hamiltonian 
systems based on the introduced concept of Dirac algebroid. The latter is a certain almost 
Dirac structure associated with the Courant algebroid TE* (Bm 'T*E* on the dual E* 
to a vector bundle t : E ^ M. If this almost Dirac structure is integrable (Dirac), 
we speak about a Dirac-Lie algebroid. The bundle E plays the role of the bundle of 
kinematic configurations (quasi- velocities) , while the bundle E* - the role of the phase 
space. This setting is totally intrinsic and does not distinguish between regular and 
singular Lagrangians. The constraints are part of the framework, so the general approach 
does not change when nonholonomic constraints are imposed, and produces the (implicit) 
Euler-Lagrange and Hamilton equations in an elegant geometric way. The scheme includes 
all important cases of Lagrangian and Hamiltonian systems, no matter if they are with or 
without constraints, autonomous or non-autonomous etc., as well as their reductions; in 
particular, constrained systems on Lie algebroids. we prove also some basic facts about 
the geometry of Dirac and Dirac-Lie algebroids. 
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1 Introduction 

The concept of Dirac structure, proposed by Dorfman [9] in the Hamiltonian framework of 
integrable evolution equations and defined in as a subbundle of the Whitney sum TA^ © 
T*A^ of the tangent and the cotangent bundle (the extended tangent oi the Pontryagin bundle) 
satisfying certain conditions, was thought-out as a common generalization of Poisson and 
presymplectic structures. It was designed also to deal with constrained systems, including 
constraints induced by degenerate Lagrangians, as was investigated by Dirac [8J, which is the 
reason for the name. 
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The need of extending the geometrical tools of the Lagrangian formalism from tangent 
bundles to Lie algebroids was caused by the fact that reductions usually move us out of the 
environment of the tangent bundles [3] (think on the reduction to so(3, M) for the rigid body). 
It is similar to the better-known situation of passing from the symplectic to the Poisson 
structures by a reduction in the Hamiltonian formalism. 

Note that the use of Lie algebroids and Lie groupoids for describing some systems of 
Analytical Mechanics was proposed by P. Libermann [25] and A. Weinstein [36], and then 
developed by many authors, for instance [U [26l [27l [28], making use of Lie algebroids in 
various aspects of Analytical Mechanics and Classical Field Theory. 

Since a Lie algebroid structure on a vector bundle t : E ^ M can be viewed as a 
linear Poisson structure H on the dual bundle ir : E* M, a properly defined 'linear' 
Dirac structure should be viewed as a generalization of the concept of Lie algebroid. Linear 
structures of different kinds on a vector bundle can be viewed, in turn, as associated with 
certain double vector bundles. The double vector bundles, introduced in \30\ [3T] (see also 
[Ml dZj) as manifolds with two 'compatible' vector bundle structures, have been successfully 
applied in \19\ l20j to geometric formalisms of Analytical Mechanics, including nonholonomic 
constraints [101 116j . To be more precise, note first that canonical examples of double vector 
bundles are: the tangent TE, and the cotangent bundle T*E of the vector bundle E. The 
double vector bundles 

T*E* ^ E 

E* 

are canonically isomorphic (cf. [24ll35j ). In particular, all arrows correspond to vector bundle 
structures and all pairs of vertical and horizontal arrows are vector bundle morphisms. Double 
vector bundles have been recently characterized [17| in a simple way as two vector bundle 
structures whose Euler vector fields commute. 

In [201 HH]) a Lie algebroid (and its generalizations) on E has been viewed as a double 
vector bundle morphism 

e:T*E^TE* (1.1) 

covering the identity on E* . This is because the linearity of a bivector field (e.g. a Poisson 
tensor) II^ on the dual bundle E* can be geometrically expressed as respecting the double 
vector bundle structures by the induced vector bundle morphism 

: TE* TE* . (1.2) 

We obtain e as the composition of the canonical isomorphism of double vector bundles 
TZr : T*E T*E* with U^, e = fi^ o 7^^. 

An application of this approach to Analytical Mechanics, in which t : E M plays the 
role of kinematic configurations, is based on some ideas of Tulczyjew and Urbahski [32 1 133 1 131]. 

Note that we can represent the morphism (II. 2p of double vector bundles by its graph Z)^ 
in the Whitney sum bundle 

TE* = TE* ®E-'T*E* . (1.3) 

The Pontryagin bundle TE* is canonically a double vector bundle: over E* and over TM(BmE, 
and the fact that e is a morphism means that is a double vector subbundle. Moreover, 
since is the graph of a Poisson tensor (in the case when £^ is a standard Lie algebroid), the 
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subbundle is a Dirac structure on E* . This immediately leads to a generalization of the 
concept of Lie algebroid: we replace the graph with any Dirac structure D on E* which is 
linear, i.e., which is a double vector subbundle of TE* . Such an object we will cal a Dirac-Lie 
algebroid. 

As was observed already in |13| , the construction of phase dynamics associated with a given 
Lagrangian does not use the fact that the bivector field is Poisson (which, on the other 
hand, induces nice properties of the dynamics), so we will use also almost Dirac structures, 
imposing no integrability assumptions. Thus, a Dirac algebroid on E will be a linear almost 
Dirac structure on E*. We introduce also afRne analogs of Dirac and Dirac-Lie algebroids. 

The main applications we propose go back again to Analytical Mechanics. To some extent, 
our concepts are similar to that of [37], where (almost) Dirac structures have been used in de- 
scription of 'implicit' Lagrangian systems. Our approach, however, we find much more general 
(we work with arbitrary vector bundles) and much simpler. This is because we obtain 'im- 
plicit Lagrangian systems' (in fact both: implicit phase dynamics and implicit Euler-Lagrange 
equations), as well as implicit Hamilton equations, just composing relations, instead of work- 
ing with somehow artificial concept of partial vector field. This generality allows us to cover 
a large variety of Lagrangian and Hamiltonian systems, including reduced systems, nonholo- 
nomic or vakonomic constraints, and time-dependent systems, with no regularity assumptions 
on Lagrangian nor Hamiltonian. 

The paper is organized as follows. In section 2 we recall basic facts concerning double 
vector bundle approach to Lie algebroids and their generalizations. Dirac algebroids, Dirac- 
Lie algebroids, and their affine counterparts are introduced in section 3, together with main 
examples. In section 4 we investigate closer the structure of Dirac algebroids, finding a 
short exact sequence of Lie algebroids associated with a Dirac-Lie algebroid and providing 
a local form of Dirac algebroids. Section 5 is devoted to inducing new Dirac algebroids by 
means of nonholonomic constraints. In section 6 we present the general schemes, based on 
Dirac algebroids, for Lagrangian and Hamiltonian formalisms. We end up with a number of 
examples in section 7 and concluding remarks in section 8. 

2 Lie algebroids as double vector bundle morphisms 

We start with recalling basic facts and introducing some notation. 

Let M be a smooth manifold and let (x"), a = 1, . . . , n, be a coordinate system in M. We 
denote with tm '■ TM — >■ M the tangent vector bundle and by ttm : T*M — > M the cotangent 
vector bundle. We have the induced (adapted) coordinate systems, {x"', x^) in TM and {x"',pb) 
in T*M. More generally, let r : — )• M be a vector bundle and let vr: E* — > M be the dual 
bundle. Let (ei, . . . ,em) be a basis of local sections of t: E ^ M and let {el, . . . ,e^) be 
the dual basis of local sections of vr: E* — > M. We have the induced coordinate systems: 
(x",y*),y* = t(e*), in E, and {x"' , , £.i = '-(cj), in E* , where the linear functions i{e) are 
given by the canonical pairing L{e){vx) = {e{x),Vx). In this way we get local coordinates 

(x^y^x^^/^■) inT^, (x^e^,x^4■) inTS*, 
{x'',y\pb,7rj) mT*E, {x-,^i,pb,^^) mT*E*. 

The cotangent bundles T*E and T*E* are examples of so called double vector bun- 
dles. They are fibred over E and E* and canonically isomorphic, with the isomorphism 
TZr'. T*E — > T*E*, being simultaneously an anti-symplectomorphism (cf. \24:\ I20j). In local 
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coordinates, TZr is given by 

7^,(x^y^pfe,7^,•) = {x%Tr,,-pb,y'). (2.1) 

This means that we can identify coordinates ttj with ^j, coordinates (p^ with , and use the 
coordinates (x", ^j) in T*E and the coordinates {x'^,£,i,Pb,y'^) in T*E* in full agreement 
with (I2.ip . According to [iTj, the double vector bundle structure is completely characterized 
by a pair of commuting Euler vector fields defining the two vector bundle structures (or by 
the pair of the corresponding families of homotheties). In local coordinates the Euler vector 
fields on T*E* read 

Vfs = Pbdp, + e.% , Vf:^ = pbdp, + y'dy. . (2.2) 

Double vector (and vector-affine) bundles will play an important role in our concepts and we 
refer to [H IZl UHl [Ml ES] for the general theory. 

It is well known that Lie algebroid structures on a vector bundle E correspond to linear 
Poisson tensors on E* . A 2-contravariant tensor 11 on E* is called linear if the corresponding 
mapping 11: T*E* TE* induced by the contraction, n(i/) = ^,^11, is a morphism of double 
vector bundles. One can equivalently say that the corresponding bracket of functions is closed 
on (fiber-wise) linear functions. The commutative diagram 




T*E 



describes a one-to-one correspondence between linear 2-contravariant tensors on E* and 
morphisms £ (covering the identity on E*) of the following double vector bundles (cf. [Ml [20]): 

T*E ^JE* (2.3) 




In local coordinates, every such e is of the form 

e{x'^,y\pb,Cj) = {x\i,,pi{x)y\c%{x)y^ik + a';{x)pa) (2.4) 
(summation convention assumed) and it corresponds to the linear tensor 

= 4(x)efc% (8) d^^ + p^(x)% d^t - a]ix)d^a ^ d^^. 

The morphism (|2.3p of double vector bundles covering the identity on E* has been called an 
algebroid in pOj. We will consider only skew algebroids, i.e., algebroids e for which the tensor 
He is skew-symmetric, i.e., is a bivector field. If is a Poisson tensor, we deal with a Lie 
algebroid. The relation to the canonical definition of Lie algebroid is given by the following 
theorem (cf. [Ill [20]). 
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Theorem 2.1. A skew algebroid structure {E,e) can be equivalently defined as a skew- 
symmetric bilinear bracket [-rje on the space Sec{E) of sections of t: E — > M, together 
with a vector bundle morphisms p: E ^ TM (called the anchor), such that 

[X,fY], = p{X){f)Y + f[X,Y], 

for f G C°°{M), X,Y G Sec{E). The bracket and the anchor are related to the bracket 
{^,ip}us = {He, dip dip) in the algebra of functions on E* , associated with the bivector field 
He, by the formulae 

ti[X,Y]e) ={i{X),i{Y)}n^, 
n*ipiX){f)) =WX),^*/K, 

where l{X) is the linear function on E* associated with the section X of E. 

3 Dirac algebroids and affine Dirac algebroids 

Let be a smooth manifold. There is a natural symmetric pairing (•[•) on the vector bundle 
TN = TN ®N T*N (called sometimes the Pontryagin bundle) given by 

(Xi + ai\X2 + a2) = ^ (ai(X2) + a2{Xi)) , 

for all sections Xi + ai, i = 1, 2, of TN = TN ©tv T*A^. Furthermore, the space Sec{TN) of 
smooth sections of TN is endowed with the Courant-Dorfman bracket, 

{Xi + ai,X2 + 02] = [Xi,X2] + Cxia2 - ix^dai , (3.1) 

where [•, •] is the Lie bracket of vector fields, Cx is the Lie derivative along the vector field X, 
and ix is the contraction (inner product) with X. An almost Dirac structure (or bundle) on 
the smooth manifold is a subbundle D of TN which is maximally isotropic with respect 
to the symmetric pairing (-j-). If additionally the space of sections of D is closed under the 
Courant-Dorfman bracket, we speak about a Dirac structure (or bundle) [3 [9]. 
Standard examples of almost Dirac structures are the graphs 

graph(n) = {Xp + ap^TpN :peN, Xp = n(ap)} , 
graph(a;) = {Xp + G TpN : p e N , Up = oj{Xp)} , 

of bivector fields 11 or 2- forms uj viewed as vector bundle morphisms, 

n : T*N^TN, Il{ap) = ia,U{p) , 
u : TN^T*N, Z{Xp) = -ixMP) ■ 

These graphs are actually Dirac structures if and only if 11 is a Poisson tensor and a; is a 
closed 2-form, respectively. 

Remark 3.1. A vector subbundle of a vector bundle over A'^ is often understood as a vector 
bundle over the whole base manifold A^. It is however clear by many reasons (see e.g. [171 
Theorem 2.3]) that we must consider also vector subbundles supported on submanifolds of N. 
Throughout this paper the term vector subbundle always means a subbundle of the original 
vector bundle supported on a submanifold Nq C A^. In this sense, our definitions of almost 
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Dirac and Dirac structure are slightly more general than those usually available in the litera- 
ture. By 'being closed' with respect to the bracket we clearly mean that the bracket of any two 
sections of TN, extending sections of D, does not depend over A'^o on the extensions chosen 
and gives a section extending a section of D. This uniquely defines a bracket on sections of 
D which is known to be a Lie algebroid bracket. 

Since the projection prj^y : TN — t- TA^ is the left anchor for the Courant-Dorfman bracket, 

i.e., 

[Xi + ai, /(X2 + a2)l = flXi + ai,X2 + + Xi(/)(X2 + aa) , (3.2) 

it is a straightforward observation that the bracket of extensions of sections of a subbundle 
D, supported on a submanifold A'^o of N, does not depend on the extensions if and only if 

prTiv(^) CTTVo- (3.3) 

Indeed, if / is on Nq, by (j3.2p Xi{f) must be on A^'o for any section Xi + ai which 
belongs to D along A'q. The condition (13. 3p we will call the first integrability condition for the 
Dirac-Lie algebroid. Under this condition the Courant-Dorfman bracket restricts to 

I-, -Id : Sec(Z)) X Sec{D) Sec{TN) . (3.4) 

Then, the second integrability condition says that [•, -Jd takes values in Sec{D): 

[■, -Id : Sec{D) x Sec{D) ^ Sec{D) C Sec((rA^)|jVo) , (3.5) 

which, according to (j3.2p and (j3.3p . is sufficient to be checked on a generating set of sections 
of D: 

{(^kiO'ijD G Sec{D) for {ai} C Sec(D) generating D. (3-6) 

By definition, an almost Dirac structure is a Dirac structure if and only if it satisfies both the 
integrability conditions, (j3.3p and (|3.6p . 

Remark 3.2. Suppose that an almost Dirac structure D satisfies the first integrability con- 
dition, i.e., the Courant-Dorfman bracket |-, -Jd of sections of D supported on A'^o is well 
defined. If we have chosen a subbundle K of TN complementary to D over A'o, we can define 
the bracket 

[•, -l^ : Sec{D) x Sec{D) Sec{D) (3.7) 

by projecting the value of [•, -Jd onto Sec{D) along K. Of course, if D is a Dirac structure, 
[•, -J^ does not depend on the choice of K and is just the Lie algebroid bracket on sections of 
D. 

In Geometric Mechanics there is often a need to use affine bundles and affine versions of 
algebroids [Ml El [HI El [El EH |22] {affgebroids in the terminology introduced in [Tn[T2]). 
We will use the following concept. 

Definition 3.1. Let A be an affine subbundle of a Lie algebroid E ^ M with the bracket 
[•, •] and the anchor p : E TM, supported on a submanifold S C M. Let V = \/{A) be its 
model vector bundle viewed as a vector subbundle of E. We call A an affine Lie subalgebroid 
in E, if the brackets of sections of A lie in Sec(v(^)), i.e., p{A) C T5 (thus the bracket of 
sections of A is well defined over S) and [cr, a'] € Sec(y) for all a, a' £ Sec{A). 
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For a more extensive treatment of brackets on affine bundles we refer to |1H [T2] (see also 
[291 mi [22]). 

To consider also affine versions of (almost) Dirac structures, we propose the following 
(compare [Tr,T2]). 

Definition 3.2. An affine almost Dirac structure on a manifold N is an affine subbundle 
D of TN, supported on a submanifold A'^o of A^, whose model vector bundle v(D) C TN 
(canonically represented by a subbundle of TN) is an almost Dirac structure on N. An affine 
almost Dirac structure is called affine Dirac structure, if the Courant-Dorfman bracket of 
sections of D makes sense (like the analogous concept for Dirac-Lie algebroids) and takes 
values in the set of sections of v(L'), i.e., ()3.3p is satisfied, so that ()3.4p is well defined and 

[•, -Id ■■ Sec{D) x Sec{D) Sec(v(i?)) C Sec{TN) . (3.8) 

The following is straightforward. 

Proposition 3.1. If D is an affine Dirac structure, then v(D) is a Dirac structure. 

Let now -F be a vector bundle over a manifold M. Since both, TF and T*F, are canonically 
double vector bundles, their Whitney sum carries a structure of canonical double vector bundle 
as well. From the general theory we easily derive the following (cf. [24t il7j ). 

Theorem 3.1. If F is a vector bundle over M, its Pontryagin bundle TF = TF (Bf T*F, 
canonically isomorphic to TF (Bf T*F* , is also canonically a double vector bundle structure 
with two compatible vector bundle structures: ti : TF — )• F and T2 : TF — )• TM (Bm F* ■ 

The core bundle ofTF, i.e., a vector bundle over M being the intersection of the kernels 
of the both projections, is in this case canonically isomorphic to T*M (Bm F ■ Moreover, the 
fibration 

(n, rs) :TF ^F ®m TM Sm F* 

is an affine bundle modeled on the pull-back core bundle, i.e., the core bundle T*M(Bm T over 
M pulled-back to F ©m TM ©m F* via the canonical projection F ®m TM ®m F* M. 

Definition 3.3. We call a submanifold Z) of a double vector bundle its double vector subbun- 
dle, if is a subbundle for each of the two vector bundle structures. 

Following the ideas of [T7], one can easily prove that this means that the two Euler vector 
fields defining the double vector bundle structure are tangent to D. One can also equivalently 
say that D is invariant with respect to both commuting families of homotheties defined by 
the two vector bundle structures (cf. [17]). 

Proposition 3.2. Let D be a double vector subbundle of a double vector bundle 

K K2 (3.9) 

Tl t[ 

Then, D inherits a double vector bundle structure with projections onto vector bundles Si = 
Ti{D), z = 1,2, where Si is a vector subbundle of Ki. 
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Proof. It is esy to see that the homothety h], being the multiphcation of vectors of the bundle 
Ti : K —?■ Ki by t G M, coincides on K2 with the homothety of the vector bundle K2 — >■ M. 
The submanifold D, being -invariant, has the base 5*2 C K2 which is /i^-invariant, thus is a 
vector subbundle oi K2 ^ M ^ Theorem 2.3]. □ 

Definition 3.4. A Dirac algebroid (resp., Dirac-Lie algebroid) structure on a vector bundle 
E is an almost Dirac (resp., Dirac) subbundle D of TE* being a double vector subbundle, 
i.e., D is not only a subbundle of ri : TE* E* but also a vector subbundle of the vector 
bundle T2:TE*^ TM Qm E. 

Remark 3.3. The above definition gives an analog of linearity of a Poisson or a presymplectic 
structure. 

We will consider also affine Dirac algebroids (Dirac affgebroids in short). 

Definition 3.5. An affine Dirac algebroid on a vector bundle E is an affine subbundle D of 
TE* whose model vector bundle v(D) C TE* (represented by vertical vectors tangent to the 
fibers of D) is a Dirac algebroid. An affine Dirac algebroid is called affine Dirac-Lie algebroid, 
if D is an affine Dirac structure, i.e., if the Courant-Dorfman bracket of sections of D is a 
section of v(D). 

According to proposition 13. H if D is an affine Dirac-Lie algebroid, then v(D) is a Dirac-Lie 
algebroid. 

Remark 3.4. We can consider as well other affine types of Dirac structures, defined on affine 
or special affine bundles, by considering vector-affine bundles of different types (see e.g. [18j). 
but we skip these considerations here in order not to multiply technical difficulties. 

In view of proposition [321 a Dirac algebroid D C TE* projects onto two vector subbundles: 
PhD = Ti{D) C E* and Vel^ = T2{D) C TM @m E, both based on a submanifold Md of M, 
giving rise to a single projection, 

= (rf , rf ) : D ^ Pho ®Mo Velfl C E* ®m (TM ®m E) , (3.10) 

which, according to theorem 13. 11 is an affine bundle modeled on the core Co of D pulled-back 
to Ph/) ©A/o Vel/), i.e., on (Ph/j ®M]j Vel/)) 'X-Mjj C-d- Note that the core Cd is a subbundle 
(supported on Md) of T*M (Bm E* - the core of the double vector bundle TE. 

The first component in Ph^i ®Mo Velo we will call the phase bundle and the second - the 
anchor relation (or the velocity bundle) of the Dirac algebroid D. The anchor relation is just 
a linear relation between vectors of E ('quasi- velocities') and vectors tangent to M ('actual 
velocities') and gives rise to the anchor map 

PD -.TM eM E DWId ^TMd (3.11) 
being the projection onto the first summand. 

To express linearity of an almost Dirac (or Dirac) subbundle of TE* in a more explicit 
way, consider adapted coordinates {x°',^i,x^,^j,pc,y'^) on TE*. The two commuting Euler 
vector fields are: 

Vi = padp^ + ijd^^ + y'^y^ + x^d^b , 
corresponding to the vector bundle structure over E* with coordinates (x,^), and 
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corresponding to vector bundle structure over TM (Bm E with coordinates {x,x,y). The 
corresponding homotheties read 

hl{x\C„x',ij,Pc,y') = {x\i^,tx\tij,tpcM''), (3.12) 
hl{x\£,i,x\ij,p,,y^) = {x\s^i,x\sij,spc,y^), (3.13) 

and a linear almost Dirac subbundle in TE* (Dirac algebroid) should be invariant with respect 
to both sets of homotheties. Note that the canonical symmetric pairing is represented by the 
quadratic function Q{x"', ^iiX^., ^j,Pc, y^) = PaX"' + y^^i which vanishes on Dirac algebroids. 

Example 3.1. The graph of any linear bivector field 

n = ^4.(x)6% A + p,^(x)% A , 

where c^j = —Cj^, is a Dirac algebroid: 

graph(n) = {(x^^^,i^eJ•,Pc,/) : = plix)y'' , ^ = 4(a:)y^a - Pj"(x)p„} . 

It is clear that Q vanishes on graph (11). This graph is a double vector subbundle, since the 
constraint functions 

x' - p\{x)y^ , ij - cf^{x)y% + p'^{x)Pa (3.14) 

are homogeneous with respect to the Euler vector fields Vi, V2- The phase bundle is here E* 
and the anchor relation is actually the graph of the vector bundle morphism p : E ^ TAf (the 
anchor map) given in local coordinates by p(x",y*) = [x"" , p\{x)y'^) . This means that skew- 
algebroids are particular examples of Dirac algebroids. The Dirac algebroids of this form, 
associated with a bivector field IT, we will call li- graph Dirac algebroids on E and denote Dn- 
The Dirac algebroid Dn is a Dirac-Lie algebroid if and only if 11 is a Poisson tensor, i.e., if 
and only if we deal with a Lie algebroid. 

Example 3.2. The graph of any linear 2-form 

^ = ^4(^)^fed^" A dx' + pi{x)d^i A dx' , 

where c^^ = — c^^, is a Dirac algebroid: 

graph(w) = {(x^e^,a;^4•,Pc,/■) : y" = pi{x)x'' , pa = c^fe(a:)6i* - pi{x)ii} ■ 

It is clear that Q vanishes on graph(a;). This graph is a double vector subbundle, since the 
constraint functions 

y' - pl{x)x- , Pa - c':,{x)Cki' + Pi{x)ii (3.15) 

are homogeneous with respect to the Euler vector fields Vi, V2. The phase bundle is here E* 
and the anchor relation is in fact the graph of the vector bundle morphism p : TM — )■ E given 
in local coordinates by p{x'^,x'^) = (x"' , pl{x)x'-'). The Dirac algebroids of this form, associated 
with a 2-form u, we will call uj-graph Dirac algebroids and denote D^j. The Dirac algebroid 
-Do) is a Dirac-Lie algebroid {presymplectic Dirac-Lie algebroid) if and only if oj is closed. 

Example 3.3. The canonical Dirac-Lie algebroid Dm = -Crm = -^wmj corresponding to the 
canonical Lie algebroid E = TM, belongs to the both above types. It is associated with the 

canonical symplectic form lom on E* = T*M and, simultaneously, to the canonical Poisson 
tensor IIm = T*M. In our local coordinates, the equations defining Dm are 

a;" = , ib = -Pb ■ 
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Example 3.4. Suppose we have a Dirac (Dirac-Lie) algebroid D on E ^ M. Let us consider 
the extension Eq = E x R as a vector bundle over Mq = M x M in the obvious way. Then, 
E^ = E* X R, TE^ = IE* X TM, and T*^;^ = T*E* x J*R. The subbundle Dq = D x Aq in 
TEq = TE* X TR, where Aq is the affine subbundle in TR defined by the constraint xq = 1 
in the natural coordinates (xo,xo,po) on TR, is an affine Dirac (Dirac-Lie) algebroid on Eq. 



4 The structure of a Dirac algebroid 

Let us start this paragraph with recalling that any section a : N ^ F oi a, vector bundle 
F ^ N (actually, of any fibration) is uniquely determined by its image cr{N) - a submanifold 
of F. We will denote this submanifold by [a]. 

Definition 4.1. Let K he a double vector bundle ()3.9p . We say that a section a : Ki — > K 
projects on the section a : M K2, if T2 projects [a] onto [a]. We will write a"^^ = a and call 
such a projectable. 

We say that a section a : Ki ^ K is suitable, if [a] is a vector subbundle of the vector 
bundle T2 : K ^ K2 ■ 

It is easy to see the following 

Theorem 4.1. Any suitable section a is projectable and \(J^^] is the image under T2a of 
the zero-section of ti. Moreover, the set of suitable sections, Suit(i^), is canonically a 
C°° [M) -module and the module morphism [T2] ■ 5 ^ is an epimorphism onto Sec(K2)- 

Suitable sections which project on the zero-section of the bundle Ki we will call 0-suitable. 
So the set Suito(i^) of 0-suitable sections is the kernel of the map [T2] : Suit(i^) — )■ Sec(-K'2)- 
A standard argument shows that that the C°°(M)-modules Suit(/C) and Suito(-fi') are the 
modules of sections of certain vector bundles over M, Su\t{K) and Suito(-ftr), respectively, but 
we will not go into details here. 

All this can be applied to the situation of the Pontryagin bundle over the vector bundle 

E*, 

TE* ^-^JMQmE (4.1) 



TM XT 



E* 

and easily explained in our standard local coordinates {x°',^,i,x'',S,j,Pc,y^)- The image of a 
section a of ri consists of points 

[x-,C^,x\x, C),i,{x, 0,Pc{x, e), y\x, 0) G TE* . 

This section is projectable if and only if the coefficients x^ and y^ depend on x only, 

x^ = x\x) , / = /(x) , 

thus a projects onto the section 

aix) = (x,x^(x),/(x)) 
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of TM (Bm E. Since being a vector subbundle means exactly being a submanifold invariant 
with respect to homotheties [T7|, a is suitable if the submanifold 

is invariant with respect to homotheties (I3.13p . i.e., 

x\x, s^) = x\x, , y''{x, = y''{x, , 

ij {X, SC) = sij (X,^), Pc{x, S^) = SPc{x, C) ■ 

As smooth homogeneous functions are linear, we get finally that x'^ and do not depend on 
^ {a is projectable) and that (,j and pc linearly depend on ^, 

ij{x,0 = i}{x)^i, Pc{x,0=pI{x%. (4.2) 

Recall that the section a is X + a, where the vector field on E* reads 

X = x''{x,C)d^,+ij{x,Od^^ 

and the 1-form a is 

a = pc{x, Odx" + /(x, Od^fc • 

Since linearity is measured by homogeneity with respect to the Euler vector field in the bundle, 
this implies immediately the following. 

Theorem 4.2. Let V be the Euler vector field in the vector bundle E* . A section X + a of 
Ti : TE* — > E* is suitable if and only if Cs/X = and C\ja = a. 

Such vector fields and 1-forms are sometimes called, with some abuse of terminology, 
linear. Hence, X + a\s suitable if and only if X and a are linear. This allows one to identify 
the bundle Suit(T'£'*) with Der(£') (Der(£;)* (g)A/ -E) with Der(£') being the bundle of quasi- 
derivations (or derivative endomorphisms or quasi- derivations) in E (see [23]). We will not 
go into details here. 

A fundamental observation is now the following. 

Theorem 4.3. // ai, i = 1,2, are suitable sections ofTE*, then |5i,52l cLnd d(cri|a2) cLf^ 
suitable. Moreover, ifa2 is additionally 0-suitable, then |ai,a2] o,nd 

|a2,5il - 2d (ai|a2) 

are 0-suitable. 

In particular, suitable sections of TE* are closed with respect to the Courant-Dorfman 
bracket and 0-suitable sections form a left-ideal inside. 

Proof. If Xi and Oi are linear, i = 1,2, then of course [Xi, X2] and CxiCt2 — 1x2^^1^ &s well as 
d(zxia2 + 1X2(^1)^ ^1^6 linear. To find the projection |cri,CT2l^^ in coordinates, let us write 

ai{x, = x\x)d,, + //(x)0% + y'{x)dii + g{{x)^jdx , 
a2{x,i) = ^^{x)d^b + //(x)Ci% + y\x)d(,i + gi{x)S,jdx . 

Then, direct calculations of the Courant-Dorfman bracket show that |cTi,52l^^ is represented 
by the tensor 
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If X and y are 0, we get 0. If x and y are 0, we get 

{y'fj + i'gl) (x)de. = 2(d(5i|52))2 • 

□ 

It is clear that having a double vector subbundle D, e.g. Dirac algebroid, we can consider 
suitable sections of D in the same manner. As the scalar products {a'i\a2) vanish for sections 
of a Dirac algebroid, out of theorem 14.31 we can easily derive the following. Let us fix a Dirac 
algebroid with an anchor relation Vel/) inducing an anchor map : Vel£) — > TM^). 

Theorem 4.4. If D is a Dirac algebroid satisfying the first-integrability condition, then the 
Courant-Dorfman bracket induces on the module of suitable sections of D a skew- symmetric 
bracket 

I, -jo ■■ Smt{D) X Suit(D) ^ Suit ((rs*)|ph^) 

such that 

for all f € C°°{Mo)- Moreover, if one of the sections is 0-suitable, the the resulted bracket is 
0-suitable. 

In the case when D is a Dirac-Lie algebroid, the Courant-Dorfman bracket is a Lie algebra 
bracket on Suit(T-E*) for which Suito(Ti^'*) is a Lie ideal and turns the bundles Suit(D) and 
Suito(-D) into Lie algebroids. Moreover, in this situation, as Suit(D)/ Suito(D) ^ Sec(VelD), 
we get a Lie algebroid bracket on the anchor bundle Vel d that gives rise to a canonical short 
exact sequence of Lie algebroids associated with the Dirac-Lie algebroid D. 

— y Suito(-D) — y Su\t{D) — ^ Vel/3 — ^ 0. 

In the case of a Lie algebroid E associated with a linear Poisson structure 11 on E* the Lie 
bracket of sections of E can be recognized inside the Lie algebroid on sections of Dji as the 
bracket of sections n(a) + a, associated with 'linear 1-forms' a, in coordinates a = y^{x)d(,i. 
The above theorem provides a generalization of this fact and, for each Dirac-Lie algebroid, 
describes the induced Lie algebroid structure on its velocity bundle. 

The next theorem characterizes the core bundle of a Dirac algebroid in terms of its anchor 
relation. 

Theorem 4.5. The core bundle Cd C T*M (Bm E* of a Dirac algebroid D C TE* is the 
annihilator subbundle Vel^j C T*M ®m E* of the anchor relation Vel/j C TM ©jvf E. 

Proof. To an element d G D that projects onto (ri,r2)((i) = {p,x,Vx) we can add any element 
Ux of the x-fiber of the core not changing the projections, so, due to isotropy, {vx, Ux) = for 
all Vx & (Vel£))a; and Cd C Vel^. The equality follows from the conditions on the rank. In 
coordinates, d is represented by 

d = x^dxa + /i% + y'd^i + gadx" 

and 

Ux = iidii, + Padx" . 
Since {d\d) = and {d + Ux\d) = 0, we have {T2{d),Ux) = 0, i.e. 

i>a + y'L = . 



□ 
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In order to describe the local form of a Dirac algebroid D, note first that, since an arbitrary 
Dirac algebroid D C TE* is the restriction to the phase bundle Ph/j C E* of a Dirac algebroid 
supported on the whole bundle E*, we can assume at the beginning for simplicity that Ph/j = 
E* . As the Pontryagin bundle TE* is, as the bundle over the projection 

(n, rs) : TE* ^ (TM (Bm E) ®m E* 

an affine bundle modeled on the pull-back bundle of the anchor bundle J*M®mE* (Theorem 
13. ip . we can write 

TE* ^{E* Qm'^M^mE) xm{T*M®mE*). (4.4) 

Note that the product Xjy^ in the above expression is not canonical, but it can be used to 
express the fact that we can add elements of T*M © E* to elements of E* © TxM © E^ 
and to serve for introducing local coordinates. Instead of the coordinates we have already 
used, it will be more convenient to introduce affine coordinates (x"^, r/^)in TxM © E^ and 
dual affine coordinates (x",Ci,Ci) i^i T*M © E*, so that {rf-,rp) represent linear coordinates 
in fibers of the anchor relation Vel/j and its (non-canonical) complementary subbundle 
TxM (B Ex = Vel/) ©m V, respectively, and the coordinates (CiXj) linear coordinates in 
the annihilators T*M (B E* = (Bm Vel^j, respectively. Note that represents the dual 
bundle Vel^. 

The points of D satisfy then rj^ = 0. Since we can add elements of the core Cd = Vel^j, 
coordinates Q are arbitrary. Therefore, there are sections ai of D associated with the canonical 
local basis ai of sections of Velu, r/* (cjj(x)) = 61 , which read 

ai(x^ei) = = 5f ,0,0' = 4'(^Kj-.o) • 

Due to isotropy, we have skew-symmetry c'--,{x) = —c'-,-{x). Now, we can add linear constraint 
Phx) in E* by introducing affine coordinates, say (x, x, ^^), such that Ph^i is expressed by x = 0, 
^ = 0. In this way we get the following 

Theorem 4.6. (local form of a Dirac algebroid) 

In the introduced local affine coordinates the Dirac algebroid D consists of points (x, x, rj, rj, C, 
for which 

5? = 0, e = 0, ^=0, Cfc = 4W0- (4.5) 
Moreover, c^^(x) = -4-(x). 

Let us note that the above constraints can be viewed as a common generalizations of 
(I3.14P and (I3.15P . The functions c^^ play the role of structure functions and rj = defines the 
anchor relation. We can write {rj, rj) as linear functions of variables (x, y) and {C, C) as linear 
functions of {p, ^) (with coefficients being functions of x) to derive constraints 

Tj{x,x,y) = 0, Ci{x,p,i) + c(^{x)r]^{x,x,y)^j = 0. (4.6) 

Example 4.1. For the Il-graph Dirac algebroid as described in example 13.11 we have 

7] = y, ^ = ±^ - pl{x)y'' , Cj = ij + p'j{x)Pa , Cb=Pb 

and the equations (j4.6p read 

x" - piix)y' = , + p]{x)Pa + c%{x)y'ik = , 

exactly as in ()3.14p . 
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5 Induced Dirac algebroids 

In this section we will show how appropriate linear (or affine) 'nonholonomic constraints' 
in the velocity bundle Vel^i give rise to new (induced) Dirac algebroids. These construction 
may be viewed as a generalization of the similar construction for the canonical Lie algebroid 
E = TM in [37]. 

Consider a Dirac algebroid D C TE* and let V he a vector subbundle of the velocity 
bundle VbId C TM®m E supported on S C Md C M. Let V = {T^)~'^iV) be the restriction 
of the vector bundle : D ^Velo to the submanifold V in the base, and let y° C T*M®E* 
be the annihilator of V. Of course, is supported on 5 as well and D Vel^j = Cd- Since 
T*M © E* is the core of TE*, we may add vectors from to vectors of the vector bundle 
Ti : TE* — > E* not changing any of two projections. In this sense, = V + 1^^^ is again a 
double vector subbundle of TE* which is no longer D, but still projects on V via T2, and on 
Ph/) via Ti. 

Theorem 5.1. The double vector subbundle in TE* is a Dirac algebroid on E. 

Proof. The subbundle is isotropic by definition, since V is isotropic as a subbundle of D, 
and is isotropic and orthogonal to V . The rank of this bundle is maximal, since first we 
loose rank by dim(Vel£)/y) and then we gain dimiV^ = dim(Vel£)/y). □ 

Definition 5.1. The Dirac algebroid we will call the Dirac algebroid induced from D by 
the subbundle V C Vel/j. 

Quite similarly, we can induce affine Dirac algebroids using an affine subbundle A of Vel o 
based on a submanifold S C M. Let V = m{A) be its model vector bundle viewed as a vector 
subbundle of Velz). Let us put A = {T^y^{A), and let C 1*M © E* be the annihilator 
of V. The vector subbundle A of T2 : D — > Vel/j is simultaneously an affine subbundle of 
Tj^ : D — >■ Ph£), thus vector-affine subbundle. Similarly as above, = A + is again a 
vector-affine subbundle of TE* which still projects on A via T2 , and on Ph/j via ti . Analogously 
to theorem 15.11 one can prove the following. 

Theorem 5.2. The vector-affine subbundle D"^ in TE* is an affine Dirac algebroid on E. 

Definition 5.2. The affine Dirac algebroid D"^ we will call the affine Dirac algebroid induced 
from D by the affine subbundle A C Vel/j. 

Example 5.1. Consider a Dirac algebroid Dji on a vector bundle t : E ^ M, associated 
with a linear bivector field 11 on i?*. Since the anchor relation Vel/jj^ is in this case the graph 
of the anchor map p : E ^ TM, subbundles V of Vel/jjj may be identified with subbundles Vq 
of E,V = {p{v)+v;v eVo}. 

It is convenient to see all this in local coordinates {x°',S,i,x^,^j,Pc,y^) in TE*. We may 
choose local coordinates in (x"") = (x°,x'^) in M, so that S is given locally by x^ = 0. Let 
us also use linear coordinates (y*) in the fibers of E, so that y = (y*) = {y'',y^) and the 
subbundle Vq is defined by the constraint y^ = 0. On TE* we have then local coordinates 
{x'^,x'',ii,pc,y'-,y^) where we have also decompositions {(,k) = {(,k,Ck) and (ii) = (^a,^l) 
associated with the decomposition (y*) = (y^y^). The double subbundle V is defined by the 
constraints (cf. example 13. ip 



V = {(x^ei,x^4■,Pc,/) : = 0, y^ = 0, = p':{x)y\ 4 = 4ki^)y'^i - Pt(^)Pa} 
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The points {x"',pb, ii) ofT*M®ME* belong to y° if and only if j;^ = and Pbp\{x)y''+iiy'' = 
for all {y''), thus .^^ = —p\{x)ph and ^/ are arbitrary. As the first condition agrees with the 
original constraints, we get the final constraints defining D^^: 

= , i'' = p\{x)y' , = c>Kix)y'Cj - pl{x)Pa , / = , (5.1) 

as adding makes S,k arbitrary. 

Let us assume now that 11 is a Poisson tensor, i.e., Dn is a Dirac-Lie algebroid. The first 
integrability condition for D^^ is now prj^, C TS*^, i.e., = (x",0)y'' = for all y\ thus 

pf(x°,0) = 0. (5.2) 

To check the second integrability condition, let us note first that D^^ is locally generated by 
the sections and R^^ = n(d^J + d^^. Since, by the assumption that 11 is Poisson, 

/th components of d (n(d^t,d^')) must vanish along E*^, i.e., 

cf,,(x°,0) = 0. (5.3) 



The vector fields d^j commute, so it remains to check whether \d^j,R£^^\ are sections of 

ia^„i?^j = [%,n(dej] (5.4) 



over E*g. But 



and, according to ()5.2I) . 

n(d6)(x°,o,Ci) = c^::(x-,o)c."% + /'(x-,o,e^)%, 

for some functions /^'. Hence 

[ag,,i?gj(x",o,e.) = %(/'')(x",o,e^)ag^, , 

so that the expression in (j5.4p is spanned over by d^^ . Thus we get that D^^ is a Dirac-Lie 
algebroid if and only if (|5.2p and (|5.3p are satisfied. This, in turn, means that Vq is a Lie 
subalgebroid in the Lie algebroid on E associated with H, so F is a Lie subalgebroid of the 
Lie algebroid Velj^ ~ the graph of the anchor map. 

Theorem 5.3. // Dn is a U-graph Dirac-Lie algebroid, then is a Dirac-Lie algebroid if 
and only if V is a Lie subalgebroid of Velz). 

Example 5.2. A particular case of the above example is the canonical Dirac-Lie algebroid 
D]\f. In this case we recover the induced Dirac structure considered in [37], i.e., the set 



Dv = {{X + a) e TT*M Sj-m T*T*M; 

X G (TTTMrHVo), yW£ (TTTMrHVo) : {a,X) =ujm{X,W) }, 

where ljm is the canonical symplectic form on T*M. Let us show that this indeed is the case. 
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According to our definition the canonical Dirac-Lie algebroid D m on the cotangent bundle 
is given by the canonical Poisson structure Hm or the canonical symplectic structure um on 
T*M, i.e., 

Dm = graph(nM) = graph(a;A,/). 

The velocity bundle Vel/j^^ C TM (Bm'^M is in this case the graph of the identity map on TM, 
the phase bundle P^Dm is the whole cotangent bundle T*M and the core C T*M©jvfT*M 
is the graph of the minus identity map on T*M. 

Any subbundle V of the velocity bundle is given by a subbundle Vq of the tangent bundle 
TM and is of the form V = {v + v; v £ Vq}. Then we get 

y = {^2") ' {V) = {X+ iOM{X) e TJ*M : TTTMiX) € Vq }. 

The anihilator consists of all pairs of covectors {if, if)) at the same point in M such 
that (/? + 'i/' G {^o)^ ■ Since the induced Dirac structure is D\j = V + , we have that 

D\j = {{X + ^) + {Cjm{X) + V-); T7rAf(X) eVo^^ + il^e {V^f }. 

The " +" sign in brackets in the above formula stands for adding an element of a core to an 
element of double vector bundle. To compare D\;[ with the Dirac structure considered in |i37j 
let us observe, that the projection of on TT*M gives the whole (Tvrjvf )~^(Vb). Adding 
elements of a core of a double vector bundle does not change projections, therefore adding 
to X produces another element Y of (T7rM)~^(Vb). Since Cjm is a double vector bundle 
isomorphism, it respects the structure of the double vector bundle. In particular, it maps the 
core of TT*M to the core of T*T*M. Both cores are isomorphic to T*M, and ujm restricted 
to the core is the identity map. We have 

CjM{y) +i^= ^m{X + ip)+^ = UJm{X) + (v9 + V) , 

so 

dX, = {X + {Cjm{X) + V); TvtmI^) G 1^0, g (Yof } • 
Evaluating u)m{X) + -0 on any W € (T7rj\,/)^^(Vo), we get that 

{CjM{X) + ij,W) = {Cjm{X),W) + {^1:J^m{W)) = {Gjm{X),W ) = um{X,W) , 

thus D\j C Dy- For dimensional reasons the inclusion is in fact equality. 

6 Lagrangian and Hamiltonian formalisms based on Dirac 
algebroids 

6.1 Implicit differential equations 

Let us start with an explanation what we will understand as implicit dynamics on a manifold 
N. 

Definition 6.1. An ordinary first-order (implicit) differential equation (implicit dynamics) 
on a manifold N will be understood as a subset T) of the tangent bundle TA^. We say that a 
smooth curve 7 : M — > (or a smooth path 7 : [tc^i] — ^ A^) satisfies the equation T> (or is a 
solution ofT>), if its tangent prolongation 7 : M — > TA^ (resp., 7 : [toi^i] — ^ TA^) takes values 
in T>. A curve (or a path) 7 in TA^ we call admissible, if it is the tangent prolongation of its 
projection 7jv on A^. 
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According to the above definition, solutions of an implicit dynamics V on a manifold N 
are projections of admissible curves 7 lying in D. Note, however, that different implicit 
differential equations may have the same set of solutions. First of all, if D is supported on 
a subset A^O; tn{T^) = Aq, only vectors from V n TA'o do matter, if solutions are concerned. 
Hence, T>' = T> r\ TAq has the same solutions as P, and T> C TAq is the first integrability 
condition. Of course, replacing V with D' may turn out to be an infinite procedure, but we 
will not discuss the integrability problems in this paper. 

All this can be generalized to ordinary implicit differential equations of arbitrary order. 
In this case we consider V as a subset of higher jet bundles, the n-th tangent bundle T"A in 
case of an equation of order n, and consider 7 as a solution when its n-th jet prolongation 
takes values in V. If we call the n-th jet prolongations admissible in T"A, then solutions of 
T> are exactly projections 7Ar to A of admissible curves (or paths) 7 in T"A lying in T>. 

Remark 6.1. The implicit differential equations described above are called by some authors 
differential relations. Let us explain that we use the most general definition, not requiring 
from V any differentiability properties, since in real life the dynamics V we encounter are 
often not submanifolds. This generality is also very convenient, as allows us to skip technical 
difficulties in the corresponding Lagrangian and Hamiltonian formalisms. Of course, what is 
a balast in defining implicit dynamics can be very useful in solving the equations, but in our 
opinion, solving could be considered case by case, while geometric formalisms of generating 
dynamics should be as general as possible. Note also that for any subset Aq of a manifold 
A the tangent prolongations TAq, T^Aq, etc., make precisely sense as subsets of TA, T^A^, 
etc. They are simply understood as families of the corresponding jets of appropriately smooth 
curves in A which take values in Aq. 

Admissibility of a path in TA^ has a natural generalization for paths 7 in an algebroid E. 
This concept plays a fundamental role in the 'integration' of Lie algebroids to Lie groupoids \2l 
and appears as natural consequence of the algebroid version of the Euler-Lagrange equations 
[13l[10]. We propose the following extension of this concept to Dirac algebroids, which reduces 
to the standard definition for Il-graph Dirac algebroids and Lie algebroids. 

Note first that given a smooth curve or path 7 with values in E we have a unique 'tangent 
prolongation' of 7 to a curve (or path) 7 : R — )■ TM (Bm E (resp., 7 : [toi^i] — ^ TM (Bm E), 
defined in obvious way by 

7(t) =7A/(t)e7(t), (6.1) 
where 7m is the projection of 7 to M, 7jv/ = r o 7. 

Definition 6.2. Let D he a Dirac algebroid on t : E ^ M and let Vel/j C TM ®m E he its 
anchor relation. We say that a curve 7 : M ^ (or a path 7 : [tc^i] — ^ E) is D -admissible., 
if its tangent prolongation -7 takes values in Vel/j, 

V t G M [7(i) = 7M(t) + 7(t) G Vel^ C TM ®m E\ . 

Remark 6.2. It is easy to see that in the case of a Il-graph Dirac algebroid, when Vel/j 
is the graph of the anchor map p : E ^ TM, a curve 7 in is D-admissible if and only 
if p{'^{t)) = ^'hiit) that coincides with the concept of admissibility for Lie algebroids. In 
particular, for the canonical Lie algebroid E = TM and the corresponding canonical Dirac 
algebroid Dm, curve 7 in TM is Djvf -admissible if and only if it is admissible, i.e., it is the 
tangent prolongation of its projection on M, 7(f) = jAiit)- 
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6.2 Phase dynamics, Hamilton, and Euler-Lagrange equations 

Our experience in working with (constrained) systems on skew-algebroids |13^ [TO] suggests us 
the following approach. Let us fix a Dirac algebroid D on a vector bundle E, 

D C T*E* ®E* TE* ~ T*E Qe* TE* . 

In generalized Lagrangian and Hamiltonian formalisms we will view D as a differential relation 

ED : T*E~\>TE* 

or 

respectively. We use the symbol ' — >' to stress that we deal with relations having domains 
in T*E or T*E* (not necessarily the whole T*E or T*E*) and with ranges being subsets of 
TE*. Note that ed = /3d°T^t and f^D is a relation over the identity on the support of D in E* 
- the phase bundle Ph/). The bundle E plays the role of the bundle of generalized velocities 
(quasi- velocities) , and its dual, E*, the role of the phase space. 

A Lagrangian function L : E ^ M and a Hamiltonian H : E* ^ M. give rise to maps 
associated with their derivatives, dL : E — )■ T*E and dH : E* — t- T*E*, respectively. The 
Lagrangian produces the phase dynamics enidL] as the image of E under the composition of 
relations = ed ° dL: 

eD[dL]=A^{E)cTE*. (6.2) 

The relation A^ we call the Tulczyjew differential of L. Similarly, when using the composition 
of relations = Pd ° dH, that projects onto the relation x% = 'Te* ° being the identity 
on a subset of E* , the Hamiltonian dynamics generated by the Hamiltonian H is defined by 

PoldH] = $g(S*) C TE* . (6.3) 

The phase dynamics e£)[dL] associated with the Lagrangian L has a Hamiltonian description, 
if there is a hamiltonian H with the same dynamics, e£)[dL] = /3/)[d//]. 

Of course, the actual phase spaces associated with L and H are projections of the phase 
dynamics on E* , Ph^ = TE*{eD[dL]) and Phg = TE*{PD[dH]). 

Since, as easily seen, the projection of the relation = sd ° dL to E ®m E* is actually 
a function, 

= TE* o Af3 = T*r o dL|v^|L , 

called the Legendre map associated with the Lagrangian L. The domain of the Legendre map 
will be denoted Vel|3 and called the Euler-Lagrange domain. It is easy to see that Ph^ is the 
image of the Legendre map 

X^ = te*oA^:Ed Vel|3 Ph^ C E* , 

and the Legendre map is the restriction of the vertical derivative d^L : E ^ E* to 

Vel|3 = {ve pr^(Veli3) : dU(t;) G Ph^} . (6.4) 

In local coordinates, 

d''L{x,y) = (^x, ^iix,y) 
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If Z) is a Il-graph Dirac algebroid, then Vel^ = E. 

The diagram picture for the corresponding Tulczyjew triple containing: T*E (the La- 
grangian side), the canonicahy isomorphic (via TZr) double vector bundle T*E* (the Hamil- 
tonian side), and TE* (the phase dynamics side) is the following (here, the arrows denote 
relations) : 




.5) 



The Euler-Lagrange equation associated with L will be viewed as an implicit dynamics on 
E. It will make sense for curves in E taking values in the Euler-Lagrange domain Vel^ C E. 

Definition 6.3. We say that a curve 7 : M — Vel|) satisfies (or is a solution of) the Euler- 
Lagrange equation, if 7 is A^-related to an admissible curve 7 in TE* (i.e., 7 is the tangent 
prolongation of its projection je* onto E*). In particular, 7 is A^-related to the curve je* 
which satisfies the phase equation. 

To describe the Euler-Lagrange equation explicitly, consider the tangent prolongation of 
the relation A^, 

TA^ = Tez) o TdL : TE — oTTE* . 

In TTE* we can distinguish holonomic vectors, i.e., vectors G T„Tii^* such that v equals 
the tangent projection of onto TE* , i.e., v = TTE*iXy). The set of holonomic vectors can 
be seen as the second tangent bundle T^E*. We define the (implicit) Euler-Lagrange dynamics 
as the subset of TE defined by the inverse image 

Sj^ = (TA^)-\T^E*)cTE. 

Theorem 6.1. If a curve j : M ^ E satisfies the Euler-Lagrange equation, then its tangent 
prolongation takes values in £j^. In particular, 7 is D -admissible. 

Proof. Let 7 be an admissible curve, 7 = 7^;*, contained in e£)[dL] and A|3-related to 7. 
Then, its tangent prolongation 7 is TA^-related to the tangent prolongation 7 of 7. But 7 is 
the 2-tangent prolongation of 7£;*, thus lies in T'^E*. □ 

Note that the converse is 'almost true'. Indeed, if 7 lies in T(A^)~^(T^£^*), we only need 
to know that we can pick up a curve in T'^E* being TA^-related to 7. This can be assured, for 
instance, by some smooth transversality assumptions. As we do not want to consider these 
questions here, let us only mention that the converse of theorem 16.21 is always true in the case 
when A^ is a map, for instance for Il-graph Dirac algebroids. 

Remark 6.3. Let us observe that in our setting the Euler-Lagrange equation is a first-order 
equation on E, in full agreement with the fact that the Hamilton equation is first-order as 
well. In the standard setting, the Euler-Lagrange equation is viewed as second-order, but for 
curves in the base M. This can be explained as follows. The solutions of the Euler-Lagrange 
equations are always D-admissible. In the case of the canonical algebroid E = TM the 
admissible curves in TM are exactly the tangent prolongations of curves in the base M, thus 
we may view the corresponding Euler-Lagrange equations as first-order equations on tangent 
prolongations, so second-order equations for curves on the base. 
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6.3 Hyperregular Lagrangians 

Let us assume that we have a hyperregular Lagrangian L : E ^ M., i.e., such a Lagrangian 
that its vertical derivative £ = d^L : E E* is a diffeomorphism. For instance, L can be 
of mechanical type, being the sum of a 'kinetic energy' (associated with a 'metric' on the 
vector bundle E) and a potential (a basic function). It is well known [T3] that in this case the 
Hamiltonian H : E* defined by 

H = {Ve{L)-L)oQ-\ (6.6) 

where V^; is the Euler vector field on the vector bundle E, defines the same Lagrangian 
submanifold in T*E* as L in T*E, when we identify canonically both bundles: 

dH{E*) = TZr{dL{E)) . 

In local coordinates, = §^{x,y) and 

H{x,0=Ci-y'{x,0-L{x,y{x,0). 

It is then easy to see that the Legendre map Xj^ is a diffeomorphism of Vel^ on Ph^, and that 
the phase dynamics associated with L and H coincide. 

Theorem 6.2. If L is a hyperregular Lagrangian, then, for any Dirac algehroid D, the phase 
dynamics e£)[d-L] coincides with the phase dynamics j^DldH] for the Hamiltonian H defined 
by \6. 6\) . In this sense, for hyperregular Lagrangians, the Lagrangian and Hamiltonian for- 
malisms are equivalent. 

6.4 Constraints 

Nonholonomic linear (or affine) constraints in our Dirac algebroid setting are understood as 
represented by vector (affine) subbundles V of the the velocity bundle Vel/j. This could look 
strange for the first sight, but it becomes quite natural, if we recall that the solution of the 
Euler-Lagrange equations are admissible curves 7 in the bundle E of quasi-velocities. Since 
there is a canonical tangent prolongation 7 of 7, with 7 lying in Vel/5, the constraint V gives 
us equations for 7 with 7 in V . The general principle is the following. 

Definition 6.4. (Nonholonomic constraints) The phase dynamics and the Euler-Lagrange 
equations for a constrained Lagrangian system on a Dirac algebroid D over a vector bundle 
T : E ^ M, and associated with the Lagrangian L : E M and the linear (affine) constrain 
bundle V C Vel^, is the dynamics associated with the same Lagrangian but on the induced 
Dirac algebroid over E. 

Another type of constraints we can consider in our setting are vakonomic constraints 
represented by a submanifold (not necessary an affine subbundle) C of E. Let us recall that 
with any submanifold C of E and any function L : C — )• M we can associate a lagrangian 
submanifold [dLc] of T*E defined by 

[dLc] = {r]eT;E:y£C and \/v G TyC (r?, v) = {dL{y),v)} . (6.7) 

We can view [dLc"] as a relation [dLc] : E — l>T*ii^. Now we can define the constrained phase 
dynamics and the Euler-Lagrange equations completely analogously to unconstrained ones, 
but replacing the relation dL{E) with [dLc]- 
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Definition 6.5. (Vakonomic constraints) The phase dynamics for a constrained Lagrangian 
system on a Dirac algebroid D over a vector bundle t : E ^ M, associated with the La- 
grangian L : E ^M. and a vakonomic constraints represented by a submanifold C oi E, is the 
dynamics represented by the subset e£)([dLc]) of TE*. We say that a curve 7 : M ^ [dLc] 
satisfies the vakonomically constrained Euler- Lagrange equation, if 7 is e£)-related to an ad- 
missible curve in TE*. 

Remark 6.4. Note first that, by definition, the phase dynamics for vakonomic constraints 
depends on the restriction of the Lagrangian L to C only. Second, we recover the old dynamics 
in the unconstrained dL{E) = [(ILe]- This can look strange at the first sight that 

we define solution of the Euler-Lagrange equations as curves in [dL^] and not in Vel^) C E, 
but when the constraints are absent there is no real difference between Velu and [dLveij,], 
since the projection tte establishes a diffeomorphism. In the presence of a constraint we have 
no longer this diffeomorphism. Of course, we could say that a curve in Vel/j satisfies the 
constrained Euler-Lagrange equation, if it is a projection of an appropriate curve in [dLc], 
but our approach seems to be more natural. It could happen that one curve is the projection 
of different curves in [dLc] that is a geometric interpretation of the presence of 'Lagrange 
multipliers'. 



7 Examples 

Example 7.1. (Mechanics on a general Dirac algebroid) 

The very general scheme of the phase or the Euler-Lagrange dynamics on a Dirac algebroid 
D C TE* can be described in local coordinates as follows. Let us choose the standard adapted 
coordinates (slightly reordered) (x, ^, x, y,p, ^) in TE* . Starting with a Lagrangian L : E —^M 
we can define the associated subset [[dL]] in TE* as consisting of points with coordinates for 
which = ^{x,y) and p = —^{x,y). Next, we intersect [[dL]] with D getting the (implicit) 
Euler-Lagrange equations defined by the following relations (in coordinates of (j4.6p ): 

x,^{x,y)] e Pho , ri{x,x,y) = 0, (7.1) 



dy 

Similarly, starting with a Hamiltonian H : E* and defining the subset [[d-fT]] by putting 
the constraints y = ^(x,^) and p = ^(x,^), we get after intersecting with D the following 
(implicit) phase dynamics 

(x,OGPhD, (x,0] =0, (7.3) 



x,^{x,i),dj (^'^'^(^'O) = 0. (7.4) 



For the canonical Dirac algebroid Dj^i we have in adapted coordinates f]^ = x"" — y"" , Ca 
+ Pa, and c^j = 0, so we get the standard Euler-Lagrange 

-df=y^ dt[dr)^''^y^ = d^^^'^y^ 
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and Hamilton 

equations. Changing the symbols y, for velocities and momenta into the standard ones, x,p, 
we end up with the traditional Euler-Lagrange and Hamilton equations. 

Example 7.2. (Pontryagin Maximum Principle for general Dirac algebroids) 

Starting with a general Dirac algebroid as above, let us impose a vakonomic constraint C C E 
parametrized by / : M x [/ — > C, with U being a manifold of 'control parameters'. In 
local coordinates {x,y) in E and u in U, the parametrization yields y = f{x,u). Note that, 
classically, for E = TM and y = x, the constraint C represents the differential equation 
X = f{x,u). 

A Lagrangian L : C — )• M may be now seen as a function L : M x [/ — > M, and [dLc] 
consists of points {x,y,p,£,) G T*E (we skip the indices) such that 

y = f{x,u), p = ^ C^ = ^- (7.5) 

ox ox Ou ou 

The above identities define a subset [[dLc]] in TE* which, similarly as above, leads to implicit 
Euler-Lagrange equations 

(x,OGPhD, ^{x,xJ{x,u))=Q, (7.6) 

{x,i^ - ^i^^y)^i) + c'i,{x)r,'{x,±J{x,u))^ = 0, (7.7) 

constrained additionally by 

^^_^ = 0. (7.8) 
du du ' 

Let us note that equations (j7.6p and (j7.7p are the same as the Hamilton equations (|7.3p and 
(j7.4p with the Hamiltonian 

H{x,u,0=^-f{x,u)-L{x,u) (7.9) 

depending on the parameter u. Moreover, the equation (|7.8p reads ^(x,m,^) = that is an 
infinitesimal form of the Pontryagin Maximum Principle (PMP): our solutions choose control 
parameters which are critical for the Hamiltonian. The whole picture is an obvious gener- 
alization of (PMP), this time for Dirac algebroids, of course in its smooth and infinitesimal 
version. 

Example 7.3. (Mechanics on skevi^ algebroids) 

Consider the Dirac algebroid Dji associated with a linear bivector field H on as described 
in example 13.11 Since in this case Djj is the graph of the map H, the relation s^jj is a map. 
Hence, A^^ = ^Du ° is also a map : E — )• TE* which in local coordinates reads 

(f)T Fit r)T \ 

x-,^{x,y),pi{x)y',4{x)y^^{x,y) + p';{x) — {x,y)j . (7.10) 

The Legendre relation A^^ is also a map which reads 



dL_ 

' dy^ 
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Let j{t) = {x{t),y{t)) be a smooth curve in E. Since 7 = o 7 is the only curve in TE* 
which is A|j,^ -related to 7, the latter satisfies the Euler-Lagrange equation if and only if 7 is 
admissible, i.e., 7 = 7. In local coordinates, 

^i-) = pti-)y'^ ^(0)(^'^) = 4(-)2/^^(-'?/) + ^|(-)^(-'^)' (7-12) 

in the full agreement with the Euler-Lagrange equation for Lie (and general skew) algebroids as 
described in [13l [lOl [261 [271 [36] . Note that we do not assume any regularity of the Lagrangian. 
As for the Hamilton equations, let us note that also in this case the relation Pd^^ is a map, 

^(x^e,,m,y^) = {x-,Cj,pi{x)y'',4{x)y'Ck - p';{x)Pa) . (7.13) 

The corresponding phase dynamics is explicit and associated with the Hamiltonian vector 
field 



XHix,0 = ( cU^)^kQ^^i^,0 - P'Six)g^ix,Oj d^^+p\{x) — {x,i)d,, , (7.14) 



I.e., 



\X)c,k 



i3={4j{x)ik^{x,0-p'j{x)-Q^Sx,i)j x' = p\{x)^^{x,i). 



In the particular case of the canonical Lie algebroid E = TM, we can take for coordinates 
y in the fiber the coordinates x'^ induced from the base. As now c^^ = (coordinate vector 
fields commute) and = 5'^ (the anchor map is the identity), we get the traditional Euler 
Lagrange equations 

dx'' .„ d ( dL \ , . , dL 



—— = x"-, — - — (x,x) = — — (x,x) , 
dt dt yOx^J ^ ^ dx''^ ' 

as a particular case. Also The Hamilton equations become completely traditional in coordi- 
nates ^ replaced by the corresponding momenta: 

Example 7.4. (Mechanics on presymplectic manifolds) 

Consider the Dirac algebroid associated with a linear 2-form uj on E* , as described in 
example 13.21 Since in this case is the graph of the map u) : JE* 'T*E* ~ T*S, the 
implicit phase dynamics associated with a Lagrangian and a Hamiltonian are inverse images 
of the images of dL and dH, respectively. In coordinates, 

/3^[dF] = |(x^enx^4■) ■■ F^a{x)x- = ||(^,0 , 4bix)^kx' - p^(:r)^, = ^(^,o} 

and 

f ■ r dL ■ dL 

eD[dL] = Ux'',Cux\ij) -.By ^i = —{x,y), pi{x)x'' = y\ c^^hix)^kx'' - pi{x)i = —{x,y) 

The implicit Euler-Lagrange equations (Euler-Lagrange relations) take the form 

Pl{x)^{x)=y\ p*^(x)^(^)(x,y) = 4(x)^(x)^(x,y)-^(x,y). (7.15) 
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Of course, for the canonical symplectic structure um = dpa A da;'* on E* = T*M we get 
the classical dynamics as above. But also in the case of a regular presymplectic form of rank 
r, 

w = ^ dpa A dx"- , 



a<r 



we get the equations for the presymplectic reduction by the characteristic distribution to the 
reduced symplectic form: the coordinates and x'^ with a > r are simply forgotten, 

d f dL\ dL 

{x,x) = —- — x,x , a<r. 



dt \dx°- J dx'^ 
Example 7.5. (Non-autonomous systems) 

Consider the affine Dirac algebroid Dq on Eq = E x M described in example 13.41 for the 
Il-graph Dirac algebroid D = Dn on E, as in example 17.31 In coordinates. 

Do = {(x^x^e„x^x^4•,po,Pc,/) : = i , = , 4- = - Piwpa} . 

For a Lagrangian L : x M — ?> M we get the Tulczyjew differential A^^ : Eq — l>Tii^Q of L 
which is the map which in coordinates reads 

such that 

= — , xO = 1, = , 4- = 4.{x)y^ — {x^,x-,y')+p>^^{x)—^{x^,x^f) . 

Identifying with the time parameter t, we get the corresponding Euler-Lagrange equations 
in the form 

^ = Pi(-)y'^ ^t{l^(t,x^,y^)^=4{x)y^^^^^ 

This is exactly the Euler-Lagrange equation on a skew algebroid for time-dependent La- 
grangians. Such equations have been obtained also as the Euler-Lagrange equations for affge- 
broids [HI O [2ll [22l [29] . For the canonical Lie algebroid E = TM, we get 

dx'* d ( dL , .A dL , 

—r- = X . — (t, X^X) = — (t, X, X) . 

dt ' dt V^x" 7 dx""^ ' 

Example 7.6. (Nonholonomic constraints) 

Consider once more the Dirac algebroid Dji associated with a linear bivector field 11 on i?*, 
as described in example 13. 1[ Consider also a nonholonomic constraint defined by a vector 
subbundle V oi E supported on a submanifold S C M. Using coordinates (x") = (x",x^) 
in M, so that 5 is given locally by x^ = 0, and linear coordinates (y*) in the fibers of E, so 
that y = (y*) = {y'',y^) and the subbundle V is defined by the constraint = 0, on TE* 
we have then local coordinates (x'*, x'', ^/,pc, y^), with decompositions {^k) = {^k,Ck) and 
(6) = {S,x,£,l) associated with the decomposition (y*) = {y\y^). The Dirac algebroid induced 
from Dn by the nonholonomic constraint V in local coordinates reads 

= {(x^C,,i;^4•,p„/) : x^ = 0, x^ = p',{x)y' , U = - P«(x)Pa , / = 0} . 
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The Tulczyjew differential A.k^v associated with a Lagrangian L : £^ — )• M is defined on V and 
associates with every {x'^,0,y'',0) £ V the set 

A^j.(x",0,y',0) = {(x^0,e^,x^4•)eTi?*:i;^ = p;'(x)y^ = —(x", 0, y\ 0) , 

f)T Fit 

L = 4A^)y'^{x-, 0, y\ 0) + — (x", 0, y', 0)} . (7.16) 

Note that the coordinates S^a of points from A^l^ (x", 0, y*", 0) are arbitrary. Curves A^^,- 

related to a curve 7(t) = (x"(t), 0, y''(t), 0) in F have thus arbitrary coordinates but the 
remaining coordinates, if the curve is admissible, satisfy the nonholonomically constrained 
Euler- Lagrange equations: 



X 



A 



dx'' 



0,y' = 0, —=p-{x^,o)y'-, (7.17) 



dt 



Note that a minimal integrability requirement is the first integrability condition for D^^, saying 
that pf{x'^,0) = 0. 

The constraint V is generalized holonomic if, independently on the Lagrangian, the above 
equations depend on the restriction of L to F only. Hence, cf^(rE",0) = and p'^{x°',0) = 0, 
so that V is generalized holonomic if and only if ^ is a subalgebroid of E. In the classical 
situation of a canonical Lie algebroid TM, the constraint V is generalized holonomic if and 
only if V is involutive, for instance V = TMq for a submanifold Mq in M. This is the 
traditional understanding of being holonomic. 

Example 7.7. (AfRne constraints) 

We can perform a similar procedure with an affine nonholonomic constraint instead of the 
linear one. Let us distinguish one variable y^ from y^ = {y^,y^^ such that the affine constraint 
^ C is defined by x"^ = , y^ = 1 , y^ = 0. The model vector bundle V = v{A) is as above 
and, as easily checked, the constrained Euler-Lagrange equations are 



X 



^ = 0, yO = l, / = 0, ^=p-(x°,0) + pr(x",0)y\ 

A(|^(x",0,ySO)) = (4Jx) + 4(x",0)2/^)|;^(x-,0,y',0)+ (7.18) 



dx" 



p^(x",0) — (x",0,y\0) 



exactly as in [10]. 

Example 7.8. (Rolling disc) 

To show how our method of Dirac algebroid works for an explicit constrained system, let us 
reconsider the case of vertical rolling disc on a plane studied in [37] . The position configuration 
space is the Lie group A/" = x S*^ x S*^ with coordinates {x^ ,x^ ,9^(p). The Lagrangian on 
JN in the adapted coordinates (x^, x^, 0, x^, x^, </j, 9) reads 

L(x^x^99,0,x^x^0,^) = Kn{{x^f + {x^f) + + ]^j2e\ (7.19) 
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The kinematic constraint due to the roUing contact without slipping on the plane is 

= R9 cos (f , x^ = R9 sin (f . 

Since the Lagrangian and the constraints are invariant with respect to translation with respect 
to x^,x'^,9, we have an obvious Lie algebroid reduction to = TA^/(M^ x 5^) = TM x 
which is a vector bundle of rank 4 over with coordinates 

{ip,ip,x^,x^,e) , 

associated with the basis of (global) sections (/i = d^, f2, fs, f^), where /2,/3,/4 come from 
the reductions of do, d^i, 8^2, respectively. The anchor p : E ^ JS^ is the projection onto 
TS^, and all the basic sections commute. The reduced Lagrangian we will denote also L, as 
it takes values exactly like in (j7.19p . 

The constraint subbundle y of -E is spanned by the sections ei = /i and €2 = f2 + R cos 93 • 
/s + Rsinip ■ /4, so we can use the basis ei = /i,e2,e3 = /3,e4 = f^, and the corresponding 
coordinates {(p,y) on E. The basis (61,62,63,64) induces the coordinate system in E* 

and adapted coordinates ((^j^jV'jO in T£'* and {f,£,,P,y) in T*E*. The constraint is now 
described by the equations = = 0, but we get non-trivial commutation relations 

[61, 62] = -Rcos (/? • 64 — i?sin if ■ 63. 

In other words, the corresponding Poisson tensor H on E* in the adapted coordinates reads 

n = i?(cos ip ■ ^4 - simp ■ 6.3)9^1 A % + A , 

and the Dirac structure induced by the constraints is 

^,1 = Ry'^^,3 sin (p — Ry^S_4 cos (p — p, ^2 = — Ry^ ^3 sin p + Ry^ cos (p } . (7.20) 
Hence, the nonholonomically constrained Euler-Lagrange equations (I7.17P take the form 

!/' = / = 0, ^ = 9', (7.21) 

dL dL dL 

Rsimp- y2^(v3,y\y2,0,0) - Rcosp- y^-^{p,y^ ,y'^ ,0,0) + — ((^, y\ y^, 0, 0) , 

di = -Rsimp ■ y^-^i^,y\y\0,0) + Rcosp ■ y^-^ip,y\y^, 0,0) 

Since 

x^ = y^ + Ry^ cos p , x^ = y^ + Ry"^ sin 99 , p = y^ , 9 = y"^ , 
the Lagrangian in coordinates {p, y) reads 

L{^,y\y\y\y^) = \m{{y^)^ + {y^)^) + 

\'h{y^? + \{i^R^ + J2)(y^)^ + mRy^iy^ cos p + y^sinp). 
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and, as show the straightforward calculations, the Euler-Lagrange equations (|7.2ip take the 
form 

/ = / = 0, ^ = y\ imR' + J,)^ = 0, Ji^ = 0. (7.22) 

Going back to the original coordinates, we get finally 

x^ = Recosip, ±^ = Resinip, <f = , 6 = 0, (7.23) 

with obvious explicit solutions. 

If the phase dynamics is concerned, in view of (j7.16p . we get that er)v[dL] = Af^v-(y) is 

parametrized by {(j),y^,y'^) as follows: 

^3 = mRy^ cos (f , U = mRy"^ sincp, = y^ , 6 = 0, 6 = 0}. 

Here 6 ^-iid 6 a-^s arbitrary, but the integrability condition allows us to describe them as 
well. Let us note that the phase space Ph^ C E* is defined by the equations 

6 = /U cos 99 • 6, 6 = sin 9? • 6, (7.24) 
where /U = . Hence, the first integrability condition gives 

6 = -^6 sin (/? • = --^66 sin (p , ^ = cos • = -^66 cos (p . 

Ji Ji 

The dynamics is Hamiltonian, since the Hamiltonian 

^^'^'^^ = W^^'^^ + 2^^^' " cosy. - RU^in^f + ^(fe)' + iU?) 

defined on E* induces the dynamics f3j-)v[dH] = e^v[dL]. The equality can be checked by 
straightforward calculations. Let us only note that, since 

. dH R^ . fR^ 2 l\ , R^ , . 

y" = — = — -6 sm V? + — COS ip-\ 43 + ^6 sm ip cos ip, 

9^3 'h \J2 rnj J2 

^ dH R^ . . fR^.2 1 A ^ 

y = — = — -6 sm ip + —43 cos sm 99 + — sm ip -\ ^ , 

044 J2 J2 \'J2 m J 

imposing the conditions = 0, = of the Dirac structure, we recover the Hamiltonian 
constraints (|7.24p . 



8 Concluding remarks 

We have introduced the concepts of Dirac and Dirac-Lie algebroid as a natural common 
generalization of a skew (resp.. Lie) algebroid and a linear presymplectic structure. Aside 
from its interesting geometrical structure, Dirac algebroids, as well as their affine counterparts 
- affine Dirac algebroids, provide a powerful geometrical tool for description of mechanical 
systems by means of generalized Lagrangian and Hamiltonian formalisms. 

The kinematic configurations (quasi- velocities) form in this framework a subset of a vector 
bundle t : E ^ M and are related to the actual velocities from TM by the so called anchor 



28 



K. Grabowska, J. Grabowski 



relation, while the phase space is a subset of the dual bundle, E* . The phase dynamics induced 
by a Lagrangian or a Hamiltonian is an implicit dynamics in the phase space described by a 
subset of the tangent bundle TE*, and the associated Euler-Lagrange equations are defined 
by an implicit dynamics in E. 

We proposed a well-described procedure of inducing a new Dirac algebroid out of a given 
one by imposing certain linear constraints in the anchor relation (velocity bundle), that on 
the Lagrangian formalism level corresponds to imposing nonholonomic constraints. Since 
imposing constraints we end up in a Dirac algebroid again, our approach does not really 
distinguish between constrained and unconstrained systems, as well as between regular and 
singular Lagrangians. Since the use of algebroids already includes reductions to the picture, 
our approach covers all main examples of mechanical systems: regular or singular, constrained 
or not, autonomous or non-autonomous etc. 

The Dirac algebroid, especially the Dirac-Lie algebroids, possess a rich and intriguing 
geometrical structure whose investigations have been started in the present paper. We are 
strongly convinced that these objects, as well as their possible generalization, will allow us 
to find a proper intrinsic framework also for field theories and other areas of mathematical 
physics. 
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